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Abstract 

A variational approach is proposed to determine some properties of the adiabatic 
Holstein-Hubbard model which describes the interactions between a static atomic 
lattice and an assembly of fermionic charge carriers. The sum of the electronic energy 
and the lattice elastic energy is proved to have minima with a polaron structure in a 
certain domain of the phase diagram. Our analytical work consists in the expansion 
of these energy minima from the zero electronic transfer limit which remarkably 
holds for a finite amplitude of the onsite Hubbard repulsion and for an unbounded 
lattice size. 
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I. INTRODUCTION 

The Holstein model involves an electron-phonon coupling which yields a polaron 
ground-state provided the amplitude of this coupling is large enough compared with the 
electronic transfer integral. This result was demonstrated by S. Aubry et al. in Ref. [@ by 
assuming that the lattice modes have a negligible kinetic energy. In this adiabatic regime, 
the description of phonons is given by the static atomic displacements. During the last 
decade, the so called small polaron of the adiabatic Holstein model has been of a great 
interest to investigate some challenging problems such as: charge density wave [pHH] and 
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high critical temperature superconductivity . These studies were based on the rigorous 
Aubry's proof which has been improved by C. Baesens and R. S. MacKay in Ref. 0. 
In the anti-integrable (AI) hmit where the transfer integral of the fermionic charges is zero 
(see [0 for a revue about the anti-integrability), it is possible to find the explicit forms 
of the Hamiltonian eigenstates and to expand them numerically with respect to the charge 
transfer. 

With no electron-electron interaction, excepted the Pauli principle, a many-electron 
problem is usually reduced to find the eigenstates of one particle Hamiltonian in which 
is introduced a chemical potential to fix the charge density (see for the case of the 
Holstein model). The ground state is thus given by the product of the one particle eigen- 
states, energies of which are smaller than the Fermi level. In some cases, the interplay of 
the Coulombian interaction with a charge screening may yield an onsite Hubbard repulsion 
which breaks this method. After extending the Holstein model with the Hubbard interac- 
tion, it is yet possible to develop a variational approach where the energy functional Fad 
is issued from the bracketing of the adiabatic Holstein Hubbard Hamiltonian. Minimizing 
this functional with respect to the displacement fields, for a zero charge transfer t shows 
that some energy minima are the products of onsite polaron, onsite bipolaron and vacuum 
site. For the weak value of t, the study of the first derivatives of Fad with respect to the 
atomic displacements indicates at least one energy minimum with a displacement field which 
is similar to a AI polaron structure. In the norm meaning, the distance between this 
minimum and the AI polaron state is bounded by a charge transfer function which goes to 
zero continuously. As a consequence, in the vicinity of the AI limit, the displacement fields 
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of the many-polaron minima are proved to evolve continuously with t which implies that 
the Fad potential has some polaron-type minima in a finite domain of the phase diagram. 



II. ADIABATIC HOLSTEIN HUBBARD MODEL 

The Holstein Hubbard Hamiltonian is written as follows: 

'H = Yj ^t^o(al«i) + gniial + a^) + vui^^rii^i - uui - T ^ C'i>C'j,<7 (2.1) 

where the atomic lattice is mapped on Z'^, i.e., i G with d = {1,2,3}, the anihilation 
operators at site i G Z'^ for both phonons and fermions are Oj and Cj, respectively. The 
corresponding creation operators are written with a dag exponent. The operator is given 
by I]o-=(T,i) ^i,crCi,a and the sum is performed over the i neighboring site indexes. For 

simplicity, the phonon contribution is only one optical phonon branch with frequency ujq. 
The electron-phonon coupling amplitude is g, the onsite Hubbard repulsion is scaled by 
t; > and the chemical potential is noted u. The displacement and momentum operators 
at site (i) are given by: 

u,= + a,) (2.2) 

Pi=i^{a\-ai) (2.3) 



parameter Eq = |^ give: 



Substituting the phonon operators in Eq.f p.l]) and dividing this equation by the energy 

hujQ 



where H is an adimensional Hamiltonian with the parameters: 

In the adiabatic hmit, the coefficient 7 is neglected which is vahd at large electron-phonon 
coupling, i.e., the amplitude g is large compared with the phonon zero point energy ?ujJq/2. 
We obtain an adiabatic Hamiltonian: 

Had = E {\< + (f - + Un,^n,^ - t E ClC,,, (2.6) 

The displacement operators {wj} are now scalar variables that are noted as a vector u G S{N) 
where S{N) = TZ'^ is a real space, dimension of which is equal to the unbounded number 
of lattice sites N. The adiabatic Hamiltonian is a sum of an electronic Hamiltonian Hf,i and 
the lattice elastic energy, i.e.. Had = J2i kuf + H^i where 



Hei = E(f - + Un.,^^,^ - t E ClC,,^ (2.7) 

The energy is now written in the variational form: 

F{u,i^) = Y.\ul+<i;\H,i\ij> (2.8) 

where 1?/^ > is a normalized electronic wave function for Ne; fermion charges. It is projected 
on the usual fermion basis, i.e., >= J^ui^ul^u > where u is a Nerniultiple site-spin index 
and \e^ >= H(jo-.)gi, C*/^. |0 >. Differentiating Eg. ( |2.8|) with respect to Ut and ip^, the 
conditions for the local extrema are: 
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H,i\ij> =E\ij> 



(2.10) 



where E is the Lagrange factor due to the ip normahzation. The Shroedinger equation 
Eq.( p.lOD imphes that ip is a Hei eigenstate with the associated eigenvalue E which is the 
electronic energy. Assuming that ip is a, H^i ground state for a given u, the total energy is 
now given by the functional Fad{u) = J2i ^uf+E^Ui) which depends only on the displacement 
variables. The ip wave function may be non-unique. 

Working with the Lqo norm in the suitable space, the F functional and its derivatives 
with respect to Ui and ip^ are continue. The adiabatic potential Fad is also continue in the 
S(N) space but its first derivatives are not necessary continue as it can be shown in the AI 
limit. 



III. THE ANTI-INTEGRABLE LIMIT 

In the anti-integrable (AI) limit, the fermionic charge transfer is zero, and thus the whole 
lattice sites are decoupled. The Hamiltonian H^i is a sum of onsite Hamiltonian Hi and the 
Hei ground state is a product of onsite eigenstate vTj with the energy Ei{ui)\ 

Hi = {-^ - fi)ni + UriiiUii (3.1) 
Hi\7Ci> =Ei{ui)\7Ci> (3.2) 



The constants (/i, U) fix which type of state is the Hi ground state for a given displacement 
Ui- Selecting a set of N^i onsite states which have the lowest energy Ei provides an electronic 
ground state ip for Hei. The site i is occupied either by a bipolaron, i.e., 2 electrons with 
opposite spin or by a polaron, i.e., only 1 electron with spin up or down or else the site i 
is not occupied. Using Eq. the optimum displacement field is such that = — 1 for a 
bipolaron, or = — 1/2 for a polaron or -Uj = if the site i is unoccupied. In the case of a 
bipolaron onsite ground state with the optimum displacement Ui = —1, the Hi eigenvalues 
are (f/ — 2/i — 1, — /i — 1/2, 0). For a polaron, the Hi eigenvalues are (f/ — 2yU — 1/2, — /x — 1/4, 0) 
and for a vacuum site {U — 2fi, —fi, 0). In order to determine whether a polaron structure 
is a local Fad minimum, it is sufficient to test the Hi ground state for the different values of 
the onsite displacement Ui G {—1, —1/2,0}. As example, for /i = —3/32 and U = 3/8, the 
Hi ground state is a bipolaron for polaron for Ui = —1/2, and it is a vacuum 

for = 0, so any displacement field consisting of an assembly of Ui G { — 1,-1/2,0} is a 
Fad minimum. On the opposite, if fi = —1/8 and U = 1, the configuration which contains 
at least a displacement Ui = —1 are not stable because the bipolaron is no longer the Hi 
ground state. 

Let us note 5i the onsite spectrum gap between the ground state energy and the first 
excited state energy. This gap is non-zero (5j > excepted for some specific values of (yU, U). 
As a consequence for nearly all (/i, U) constants, the onsite ground state is not degenerate 
excepting the spin degeneracy which occures for the polaron. 
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IV. EXPANSION OF POLARON STRUCTURES 



As soon as the transfer integral t is non zero, one may guess that the AI polaron states 
should still be Fad minima, at least for a certain range of parameter. Here we propose a 
proof that confirmes this guess. Our demonstration is based on the Fad gradient study in 
the space of the displacement configurations i5(N). The Fad gradient is given by: 

<pi = Ui^ ^ (4.1) 

Let us introduce the following operators: 

P2,i = Ui^il-riii) 
P3,i = nuil-rii^) 

P4,i = l-riii- nil + fii^riii (4.2) 

They verify J2a Pa,i I >=| > for all ip state. We choose to note Pa,i \ i' >= Xa,i \ i'a,! > 
where the state ipa,i is normalized and Xa,i is a real positive scalar. The gap 6i is assumed 
to have a lower bound 6, i.e., 6i > 6 > which is valid for nearly all {fi, U) constants. Then 
only one Xa,i is non-zero at the AI limit and it is equal to one, the corresponding index a is 
noted gi. If gi is equal to either 1 or 4 then we note Xg. = Xa=g,,i- If gi is equal to either 2 
or 3, the site i is occupied by a polaron with a spin up or down. In such a case, the onsite 
ground state is spin-degenerate and we note Xg^ = \Jx\^^ + xl^. As soon as the electronic 
transfer is non-zero, Xg^ variates with t and it can be proved that (see Lemma [V^) : 

7 



1 



where ric is the number of the nearest neighbors, Ug is the number of distinct Hi eigen values. 
The inequahty Eq. ( [4 .31) imphes the continuity of Xg^i with respect to t in the vicinity of the 
AI hmit. 

The potential Fad is now derived with respect to uf. 

dFad ^ . .11. 



Ui 



Q + < V'.JylV'^,, > (4.4) 

As Pa,iP^,i = ^a^pPa^i wherc 5a^j3 is the Kroeneker symbol and as rii commutes with the P^^i 
operators: 

0i = Mi + ^ Xl^li < > (4-5) 



Let note na,i =< V'a,«l^i|V'«,j > and let write Ui = Ui{0) + pi where Mj(0) is the onsite 
displacement at t = 0. The Eg. ( [4 .51) is rewritten as follows 



0. = A-^(l-4) + ^E<.^«. (4.6) 

where we used Eq. ( p.9|) in the AI limit to find Mj(0) = The sum in the right hand 

side of Eq. ( [4. 61 ) is performed over the indexes a ^ gi for Qi = {1,4} and over the indexes 
a = {1,4} for (7. = {2,3}. 

The scalar product "p^. is now detailed: 



Focussing on the terms of the right hand side sum, some simple arguments give the following 
inequalities: 



ft 

Pi.{l-xl)^<sup,{l-xl)\p,\ (4.8) 

P^■Y.<^'^>-^^P^('^-0\P^\ (4-9) 

So each term of the sum in the Eq. ( [4.7| ) verifies 

+ p,.(-(l - + E <.^) > pI - 2|P.M«J>.(1 - 4) (4-10) 

which is positive if \pi\ > 2supi{l — x^.) and with Eq. ( ^I3| ) it is equivalent to \pi\ > Rai 
with writing: 

- (5-tn,K-l))^ ^'-^^^ 
Defining the subset -B(-R) such as u G if \ui{0) — Ui\ < R for all i index, the boundary 

of B{R) is noted B{R). For any m G > -Ra/), the product ('u(O) — u).(f){u) is positive 
and thus there is at least one displacement configuration Umm{t) G B{Rai) which is a local 
minimum of Fad- As a consequence. The potential Fad has a minimum in the u{0) vicinity 
at most at a Rai distance in the L^o norm meaning, i.e., supi\umin,i(t) — Ui{0)\ < Rai- 

V. CONCLUSION 



The present proof holds for any displacement field which is a minimum of the adiabatic 
potential Fad in the AI limit and such as the gap 5 7^ 0. For nearly all (/i, U) values, 
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the displacements of the polaronic minima variate continuously with respect to the fermion 
transfer t in the vicinity of the AI limit where t=0. So it is about the total energy of these 
minima because of the Fad continuity with respect to the displacements. Consequently, the 
adiabatic potential has some minima with a many-polaron structure in a finite region of the 
phase diagram. Nevertheless, the absolute minimum of Fad-, i-e., the ground state of Holstein- 
Hubbard model can not yet be determined for all the parameters. To that purpose, a 
numerical investigation might be required but no idea emerge to tackle the case of a non zero 
Hubbard coupling with many charge carriers, excepted to compute a meanfield as proposed in 
Ref. 0] or to study a model with a small electron number. The latter possibility is presented 
in [^,0 where the phase diagram is calculated for 2 electrons. For a two-dimensional atomic 
lattice, a critical point where coexist 3 different types of bipolaron was found far from any 
trivial limit. In this specific region, because of the bipolaron degeneracy, the bipolaron 
tunneling (or equivalently the inverse of the ground state effective mass) is very sensitive to 
the quantum fluctuations which are yielded by a non- zero 7 (Eq. ( |2.5| )) Ref. Around the 
critical point, both the bipolaron ground state mobility and its binding energy reach lOO-ft' 
with realistic input parameters (Eq. (|2.5|)). This result allowed some conjectures about the 
mechanism which yields the high critical temperature superconductivity of cuprates. 

The present study can be extended straighforwardly for a non harmonic phonon potential 
as for an atomic lattice embedded in an external magnetic field. However an extension to 



a different electron-phonon coupling such as the SSH model [|TT| might have non trivial 
anti-integrable limit which makes our arguments much less efficient. 
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VI. LEMMA 



We choose the ground state energy as the energy reference. We shall assum the non- 
degeneracy of this ground state so only one projector Pg. is such as HiPg.\ip >= for all ip. 
As \ip >— J2aPa,i\'ip > we havc \ip >— J2a^a,i\'ipa,i > where is introduced the normalized 
states ijja,i and Xa,i\ilJa,i >— Pa,i\'^ > which implies J2a \^a,if — 1- It is possible to choose 
the ipa,i such as the Xa,i are real positive for all a. One notes Ki = J^a Pa.iHeiPa,i and 
i^T = Hgi — Hi + t J2i:j,cr Ci^a^j,a- + C'^vtQ.o- whcrc (i : j) are the i neighboring sites. As 
Pa,i[Ci^a^j,(T + Cj'aCi^a]Pa,i — and as and Hi commute with the projectors Pa,i: 

Ki^H, + Hi (6.1) 

Let $0 be the Ki ground state with energy Eq. Using Pa,iP/3,i = Sa,f3Pa,i where 5^,/? is the 
Kroeneker symbol, then XjPc«,i|*o >= Pa,iKi\^o >= -EoPa,i|*o >■ 

As Hj and Hi are decoupled, the state 0o is a product of the Hj ground state and the 
Hi ground state. It follows that -PgJ$o >= l^o > and Pa^gJ^o and it is now easy to 
estabhsh that 

<$o|^/t|$o> =Eo (6.2) 
<^o\Hei\%> ^ < ^o\Ki\% >^ Eo (6.3) 
<^,Ji^#3, > >Eo (6.4) 

From the identity P^^ = Pa,i, one deduces that for all \ijja,i >= PaAi^ Pa,i\'^a,i >= 
> and < ipa,i\Ki\ipa,i >=< i'a,iWei\i'a,i >■ K (^i is the first excited state energy of Hi. 
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< iJaj^g„i\Hi\lpaj^g„i > > Si (6.5) 

< ^a^a,,i\Ki\lpa^g,,i > = < Ipaj^g^AHillpa^g^i > + < '^a^g,,i\Hi\'^a^g,,i > (6.6) 

< '^aJ^g„i\K^\i^o,J^g^,i > > < i'a^g.AHiH^a^g.^i > +Si > Eq + Si (6.7) 



The previous results Eq. (|6.7|) is necessary valid if ijj is the Hf.i ground state: 

< >> -^0 + Si (6.8) 

Multiplying H^i by the identity {J2a Pa,i = Id) gives: 



H,i = K,+ J2 Pa,iHelPa',i (6.9) 



and bracketing by ip: 



< IplHelllp >= Xl^li < > + X! ^o^,i^P,i < i^a,i\H el\ll) p^i > (6.10) 

As ?7,j o- and rij o-/ commute with each other for all suffix i, j, a and cr', and as Pa,iPi3yta,i = 0: 

^a,i^l3,i < i^a,i\Hel\tpp,i >= J2 ^»,iXf3,i < i)a,i\^\lpf3,i > (6.11) 



One deduces a simplification of the Eq. (|6.10|) : 



< 1p\Hel\lp >= "^xli < 1pa,i\Ki\i^a,i > ^a,iXp,i < V'a.i | ^ | ^/'/3,j > (6.12) 

and combining the inequalities ( |6.4|j6.8| ): 

< tp\Hei\llj >> Eo + 5iJ2 - ^ J2 < ^a,i\^\^f3,i > (6.13) 
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As we choose to map the atomic lattice on Z'^ ^ the number of first neighboring sites 
is = 2d. With Eq. (|6.3|), one now writes the set of following equations where E =< 

< ^o\Hel\^o > >E 
Eq>E -> Eo> E > Eq + 8iY^ -tY^ Xa,iXp^i < i:a,i\^\i'l3,i > 

< ?/;„,i|A|V'/3,i > = < '^taCj^a + C+^C^,^ | ^/-^^i >< 2n^ 

i:j,cr 

Y ^ ''^ Xa,iXp^i < ljja,i\/^\'lljfi,i > 

(1 - xl)6i < n^t[CY Xa,if - Y 

a a 
a 

(1 - xl)6, < n^tHY ^o., + xgf - 1] (6.14) 
The Cauchy-Schwartz inequality applied to the sum Y^a^g Xa,i gives: 

a¥=g V "^3' V 

E -^o.* — V^\/ "~ 

where = 3 is the maximum number of distinct Hi eigenvalues. With Eq. ( |6.14| ), it follows 
that 

(1 - xl)S, < net[(n|(l - xlf^ + xgf - 1] (6.15) 
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and it is now easy to obtain: 

1 

2 X 1 2nrns t 



l-<.)5<- i-^ (6.16) 



The latest result holds for the case of a non-degenerate Hi ground state. Regarding the 
case of a polaron at site i, the atomic orbital is occupied by 1 electron with either spin up 
or spin down. The onsite ground state is spin-degenerate but the same arguments as for the 
non-degenerate case can be used to establish the following inequality: 

Let write x^. = x| ^ + Xg ^ such as this inequality is now written: 



(1 - J5 < % L.— (6.18) 

^ '''' 5,-n,t{n,-2) ^ ' 

If we now assum that there is a lower bound 5 for the onsite gap (5j > 5 > 0, then the 

combination of the Eqs. ( |6.16 , 5.18 ) gives for all site i : 

1 

, o V 1 2TlrTl t , , 

{l-xl)^<- 6.19 
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